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ABSTRACT 

Monte  Carlo  calculations  are  used  to  study  theoretically  hot  electron 
transport  through  a thin  insulating  film  subjected  to  a high,  uniform 
electric  field  F.  A constant  energy-independent  mean  free  path  A leads  to 
a stable,  steady-state  energy  distribution  for  the  electrons,  character- 
ized by  an  average  steady-state  energy  E . E depends  on  X,  F and 
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Cpk»  the  optical  phonon  energy  associated  with  scattering  of  the  hot  elec- 
trons by  the  lattice.  An  energy-dependent  mean  free  path  A(E)  which  in- 
creases with  increasing  electron  energy  can  lead  to  energetic  runaway  of 
either  a relatively  small  nunfter  of  electrons  in  the  distribution  (quasi- 
stability, or  a bi-modal  distribution)  or  the  entire  distribution  (in- 


stability) . A graphical  method  has  been  developed  to  gain  insight  into 


this  problem.  The  method  rests  on  an  analysis  of  the  interactions  of  two 


curves  plotted  in  the  E,  A plane,  one  curve  being  the  plot  of  E 

ov  )8o 

where  A is  the  mean-free-path  parameter  taken  to  be  constant  (energy 


independent)  for  each  calculation  of  E , the  other  curve  being  the 
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plot  of  the  actual  functional  dependence  of  A on  E,  A(E).  The  method  is 


illustrated  by  several  prototypical  examples 


I . Introduction 


- 1 - 

In  a previous  paper,1  the  authors  studied  theoretically,  using  Monte 
Carlo  techniques,  hot-electron  transport  through  thin  insulating  films  in 
the  presence  of  high,  applied  electric  fields,  on  the  order  of  a aegavolt 
per  centimeter.  This  study  has  particular  relevance  to  contemporary  solid- 
state  electronics,  where  thin  insulating  layers,  usually  on  a silicon  sub- 
strate, are  frequently  subjected  to  electric  fields  near  the  threshold  for 
electrical  breakdown.  The  previous  paper,  hereafter  referred  to  as  I, 
studied  transport  characterized  by  a constant,  energy-independent  mean  free 
path  (mfp)  between  collisions  with  the  lattice.  The  major  result  found  is 
that  the  constant-mfp  constraint  always  leads  to  a stable,  steady-state 
energy  distribution  no  matter  how  anisotropic  the  microscopic  scattering 
process  may  be,  short  of  pure  forward  scattering.  Further,  the  average 
energy  of  the  steady-state  distribution  is  given  by  a simple  scaling  law, 
namely,  Eq.  (1)  of  the  present  paper. 

In  the  present  paper  we  study  hot-electron  transport  characterized  by 
an  energy-dependent  mfp,  namely,  one  which  increases  with  increasing 
electron  energy.  This  latter  property  is  relevant  to  the  electrical  break- 
down problem,  since  it  favors  runaway,  in  energy,  of  either  a small  number 
2 

of  electrons  or,  in  some  cases,  of  the  entire  energy  distribution.  Both 
of  these  situations  are  examined  in  the  present  paper.  In  particular,  we 
exhibit  a powerful  graphical  technique,  utilizing  the  previous  results  in 
I for  an  energy-independent  mfp,  to  analyze  the  stability  of  the  overall 
energy  distribution  where  the  mfp  is  energy-dependent,  and  further  to  obtain 
the  average  steady-state  energy  of  the  distribution  when  the  distribution 


is  stable. 
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II.  Outline  of  the  Method  and  the  Results  of  I 

High  field,  hot-electron  transport  Is  studied  In  an  aaorphous  film 
of  thickness  d subjected  to  a uniform  electric  field  F.  The  electron  Is 
Injected  Into  the  film  at  plane  z-0  In  a random  direction  and  with  an 
initial  energy  EQ.  Under  the  influence  of  the  field  It  accelerates  in  the 
positive  z-dlrectlon,  following  a Newtonian  trajectory  up  to  the  point  of 
an  inelastic  collision  with  the  lattice.  The  length  of  this  trajectory  Is 
determined  randomly  according  to  the  functional  energy  dependence  of  the 
mfp . The  inelastic  collision  process  creates  an  optical  phonon  of  energy 
e^,  and  the  scattered  electron,  therefore,  emerges  from  the  collision 
with  an  amount  eph  subtracted  from  its  energy.  The  new  direction  of  motion 
of  the  electron  is  chosen  randomly  according  to  the  functional  form  of  the 
angular  scattering  process.  He  have  taken  the  angular  scattering  to  be  de- 
pendent only  on  the  polar  angle  9 with  respect  to  the  direction  of  travel 
just  prior  to  the  collision.  For  simplicity,  we  used  a (1  + cos 0) -dependence 
to  study  dominant  forward  scattering  and  a (1  - cosO) -dependence  to  study 
dominant  backward  scattering;  isotropic  scattering  has  no  6 - dependence. 
With  its  post-collision  energy  and  direction  specified,  the  electron  again 
follows  a Newtonian  trajectory  up  to  its  nest  collision  with  the  lattice, 
and  so  on.  In  this  manner,  the  electron's  path  through  the  film  is  followed 
until  its  emergence  from  the  film  at  plane  z-d,  at  which  point  the  computer 
stores  the  emerging  electron  energy.  Another  electron  is  then  injected  at 
z-0  with  energy  E0,  and  its  path  is  followed  in  a similar  manner.  The 
process  is  repeated  until  a representative  histogram  of  the  emerging  energy 
distribution  is  obtained.  A minimum  of  200  electrons  was  used  to  obtain 
the  energy  histograms. 


r 


r 
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In  the  situation  studied  in  I,  the  film  thickness  range  was  on  the 
order  of  one  hundred  to  several  thousand  Angstroms;  Eq  was  generally  in 
the  range  0.01  to  16  eV;  F was  of  the  order  of  10^V/cm;  and  was 
either  0.05  or  0.1  eV— -all  values  appropriate  to  the  study  of  hot  electron 
transport  in  thin  films  of  A^O^  and 


The  major  result  of  I is  as  follows:  For  a constant  mfp  the  distri- 

bution-in-enecgy  for  the  injected  electrons  reaches  a steady-state  value, 
independent  of  the  initial  Injection  energy,  for  which  the  average  steady- 
state  energy,  E , is  given  by  the  scaling  law: 
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where  k is  a constant  determined  by  the  details  of  the  scattering  process: 
the  greater  the  probability  of  forward  scattering,  the  larger  the  value 
of  k.  For  dominant  forward  scattering  characterized  by  a Cl  + cosd)  - 
dependence,  k ■ 1.27;  for  isotropic  scattering,  k “ 0.87;  for  dominant 
backward  scattering  characterized  by  a (1  - cosO)  - dependence,  k - 0.66. 

A graph  of  k vs.  an  anisotropy  parameter  f,  defined  for  a simple  one- 
dimensional  random-walk  simulation,  is  presented  in  Fig.  13  of  I.  The 
'development'  distance  D,  over  which  the  electron  energy  distribution  ex- 
ponentially settles  down  to  its  steady-state  is  also  given  by  a scaling 
law:  D-k'  (FA  For  a film  thickness  d exceeding  a few  D the  emerging 

energy  profile  is  Independent  of  d. 

We  found  in  I that,  under  the  constraint  of  constant,  energy- inde- 
pendent mfp,  any  degree  of  forward  scattering,  short  of  exact  forward 
scattering,  leads  to  a finite,  stable,  steady-state  energy  distribution. 
Physically,  this  is  understood  as  follows:  For  a given  degree  of  forward 

scattering,  particles  with  energies  less  than  E are  quickly  directed  along 
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the  positive  s-direction  by  the  applied  field,  F.  This  tendency  to  travel 
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with  the  field  implies  that  particles  with  energies  less  than  E are 
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warmed  up  by  the  field.  On  the  other  hand,  particles  with  energies  greater 
than  E are  not  significantly  influenced  by  F and  travel  in  nearly 

8V  f 88 

rectilinear  trajectories.  So  long  as  there  exists  some  chance  that  an 
energetic  particle  be  back-scattered  along  the  negative  z-direction,  such 
a back-scattered  particle  will  continue  traveling  backwards,  losing  energy 
as  it  travels,  against  the  field.  In  the  case  of  strong  forward  scattering, 
this  backward  travel  is  sustained,  for  the  particle  is  now  most  likely  to 
continue  being  scattered  along  its  present  course,  which  is  the  negative 
z-directlon.  This  cooling  effect  continues  until  the  particle  is  turned 
around  by  the  field,  or  by  angular  scattering,  into  the  positive  z-direction 
Thus,  an  energetic  balance  is  struck  between  the  warming  of  low-energy 
electrons  by  the  field  and  the  cooling  of  energetic  electrons  by  back- 
scattering  into  trajectories  directed  against  the  field.  This  balance 
produces  a finite,  steady-state  distribution. 


III.  Energy-Dependent  Mean  Free  Path;  Unstable  Energy  Distributions  and 
Runaway  Electrons 

If  the  mfp  is  energy-dependent,  X ■ A(E),  and,  further,  if  X(E)  in- 
creases with  increasing  energy,  the  argument  on  stability  presented  at  the 
end  of  Section  II  breaks  down.  The  hotter  a hot  electron  gets,  moving, 
roughly,  'with'  the  field,  the  less  likely  it  is  to  make  any  collision 
with  the  lattice.  Further,  even  if  it  does  make  a large-angle  scattering, 
turning  it  into  a trajectory  'against'  the  field,  as  it  subsequently  under- 
goes cooling,  the  probability  of  its  making  another  large-angle  collision 
which  will  return  it  to  a motion  'with*  the  field  keeps  increasing.  The 
net  effect  favors  higher  energies  for  the  electrons.  Whether  the  dynamic 
interplay  of  these  effects  leads  to  full  stability,  as  in  I,  or  to  total 
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instability  of  the  distribution,  with  almost  all  electrons  increasing 
their  energy  indefinitely  with  increasing  film  thickness,  or  to  a quasi- 
stability with  most  electrons  settling  down  within  a finite  steady-state 
distribution  and  a few  running  away  in  energy  is  the  main  subject  of 
this  paper. 

We  have  studied  the  problem  of  stability,  quasi-stability  and  in- 
stability using  a graphical  technique  in  combination  with  the  results 
previously  obtained  in  I.  The  quantities  that  are  plotted  graphically 
are  E vs.  X as  given  by  (1),  where  E is  identified  with  E calculated 
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as  if  X were  constant,  and  X " X(E),  the  given  functional  dependence  of 
X on  E.  [In  plotting  the  constant  mfp  scaling  law  (1),  the  quantities  F 
and  must  be  specified,  and  the  numerical  value  of  k depends  on  the 

detailed  form  of  the  angular  scattering  law.]  The  properties  of  stability, 
quasi-stability  and  instability  are  analyzed  in  terms  of  the  intersections 
of  the  two  quantities  plotted.  Two  simple  functional  forms  of  energy  de- 
pendence of  the  mfp  were  chosen  to  illustrate  the  technique:  X ■ XQ  + 

1/2 


and  X - X + 3(E/e  .),  where  X . a and  B are  adjustable  con- 
pn  o pn  ’ o 


c,(E/tph) 

stants  with  the  dimension  of  length. 

A.  Square-root  dependence  of  X on  E 

Here  we  take  the  functional  dependence  of  X(E)  as 


1/2 

X - XQ  + a(E/eph) 

In  Fig.  1,  Eqs.  (1)  and  (2)  are  plotted  in  the  E-X  plane.  Equation 
(1)  is  a parabola  in  E vs.  X,  here  plotted  for  two  different  values  for 
eph,  namely  0.1  and  0.05  ev,  and  labelled  1 and  II,  respectively.  In 
each  case  there  is  one,  and  only  one,  intersection  between  the  X(E) 
parabola  (2)  and  the  scaling-law  parabola  (1) , given  analytically  by  the 


'co-ordinates’ : 
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Colnciding  with,  or  extremely  close  to,  these  points  of  intersection  (3) 
in  Fig.  1 are  large,  open  circles.  These  circles  represent  direct, 
independent  Monte-Carlo  calculations  of  E __  for  the  identical,  corres- 
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ponding  problem.  Thus,  a single  intersection  of  A(E)  and  E (A) 
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establishes  the  stability  of  the  distribution  and  also  yields  the  average 
energy.  Note  that,  as  with  the  A * constant  case  of  I,  there  is  not  a 
single  'runaway'  electron.  That  the  distribution  settles  down  into  a 
steady  state,  independent  of  film  thickness,  is  seen  in  Fig.  2 for  four 
different  cases,  labelled  I through  IV.  For  c "2  I the  independence  of 
E from  the  initial  energy  E is  also  explicitly  exhibited. 

SV p 88  u 

Note  that  in  the  E-A  plane  of  Fig.  1 the  constant  mfp  situation  is 
simply  the  vertical  line:  A ■ constant.  In  effect,  the  square-root  de- 

pendence of  A on  E,  given  by  (2),  is  just  too  weak  a dependence  to  give 
a result  qualitatively  different  from  the  A - constant  case  of  I. 

That  a single  intersection  of  the  scaling- law  (1)  and  the  A - A(E) 
curve  establishes  a stable  distribution  and  simultaneously  gives  E 
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is  undoubtedly  true  under  more  general  circumstances  than  studied  above. 
In  order  to  explore  this,  we  studied  a model  problem  in  which  the  angular 
scattering  is  also  energy-dependent,  namely,  one  corresponding  to  polar- 
optical  phonon  scattering  as  given  by  Eqs.  (20) -(23)  of  I,  but  with  the 
mean  free  path  having  the  weak  energy-dependence  (2).  The  graphical 
solution  exhibited  in  Fig.  3 for  two  A(E)  cases  again  yields  a single 
point  of  intersection  for  each  case.  Note  that  for  polar-optical  phonon 
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0 50  .100 

MEAN  FREE  PATH,  X (A) 

Fig.  1:  Graphical  determination  of  average  steady-state  energy  when 

the  mfp  is  a function  of  energy.  Curves  I and  II  are  E * 

av,ss 

k (FX)2/eph  with  eph  - 0.1  and  0.05,  respectively.  F - 106  V/cm 

and  k - 0.87  (isotropic  scattering).  Curves  A through  E are 

1/2 

plots  of  A - Ao  + a <E/eph)  . (A)  Aq  « 25  X,  a - 2 X; 

(®)  A - 40  X,  a - 1 X;  (C)  A - 40  X,  a - 2X;  (D)  A - 40  X, 
a - 4 X;  (E)  Aq  ■ 50  X,  « ■ 2X.  The  open  circles  represent 
Eflv , as  comPute<*  through  a Monte  Carlo  simulation  using  the 
parameters  for  scattering  and  for  A(E)  given  by  the  two  curves 
that  intersect  at  E . 
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Fig.  2:  The  development  of  Eflv  gg  with  film  thickness,  d,  for  X - X 

1/2 

a (E/Cp^)  , using  some  of  the  parameters  of  the  curves  in 

Figure  1.  F ■ 10**  V/cm  and  k * 0.87  for  all  the  curves. 

I (a):  X - 25  X,  a * 2 X,  e . ■ 0.1  eV,  and  E - 0.01  eV: 
o pn  o 

1 (b):  same  as  1(a),  except  Eq  ■ 10  eV;  II:  Xq  - 25  X, 

a - 2 X,  eph  - 0.05  eV,  Eq  - 0.01  eV;  III:  XQ  - 40  X,  a « 

Eph  - 0.05  eV,  and  Eq  - 0.01  eV;  IV:  XQ  - 40  X,  a « 2 X, 

e . ■ 0.05  eV,  and  E ■ 0.01  eV. 
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Fig.  3:  Graphical  determination  of  E when  E * k (FA)8/3/e  3 
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and  A(E)  - Aq  + a^E/e^)1/2.  k = 1.64,  F = 106  V/cm,  and 
eph  ■ 0.1  eV.  The  form  of  Eay  gs  is  generated  using  the  polar- 
optical  phonon  scattering  formula. 
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scattering,  as  ■ 1.64  (FA)8^3eptl""5^3.  The  Monte-Carlo  computation 

for  E yields,  for  each  case,  the  steady-state  energy  represented 

&V  ) ss 

by  the  heavy  dot.  It  is  seen  that  the  agreement  is  excellent. 

Qualitatively  different,  and  more  interesting,  possibilities  de- 


velop if  the  functional  dependence  of  A on  E is  stronger.  This  is 
illustrated  by  the  next  case. 

B.  Linear  dependence  of  A on  E 

Here  we  take  the  functional  dependence  of  A on  E as 


A - Aq+  0 (E/eph> • 


The  graphical  method  consists  of  simultaneously  plotting  Eqs.  (1)  and 
(4),  as  in  Fig.  4.  It  is  seen  from  this  figure  that  there  may  not  be 
any  intersections  of  the  two  equations,  or  there  may  be  two  inter- 
sections. Also  there  may  be  a unique,  single  intersection  which  is 
really  a degeneracy,  or  coincidence  of  the  two  intersections.  Alge- 
braically, these  situations  are  represented  as  follows: 


Single  (degenerate)  intersection:  6 ■ B 


-ia  Khl 

4k[FAoJ  * 


A - 2A  , E - 4k 
cr  o cr 


(FA  ) A e . 
o o ph 


where  (A  , E ) are  the  co-ordinates  of  the  single  intersection, 
cr  cr 


No  intersection: 


Two  intersections: 


Y - g6-  > 1 • 

Pcr 


-^-<1  • 


I 


Fig.  4 


Display  of  points  of  intersection  between  E - 0.87  (FA)2/e 
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and  A - Aq  + 0 (E/eph)»  wher*  E ■ 106  V/cm  and  epli  ■ 0.05  eV. 
Note  that  for  some  values  of  A_  and  8 there  are  no  intersections 
which  Implies  a runaway  condition. 


4 


where  (A+,  E+),  (A_,  F._)  are  the  co-ordinates  of  the  two  intersections. 

In  Fig.  4 the  parameter  values  are  F * 10^  V/cm,  c . ■ 0.1  eV  and 

ph 

k ■ 0.87.  For  XQ  ■ 25  X,  Eq.  (5)  gives  8cr  * 1.15  X for  the  single 

intersection  point.  Thus,  for  8 > 1.15  X there  are  no  intersections 

and  for  8 < 1.15  X there  are  two  Intersections.  Indeed,  it  is  seen  that 

for  8*2  and  1.2  X there  are  no  intersections,  and  for  8 ■ 1.0,  0.8  and 

0.5  X there  are  two  Intersections  in  each  case.  For  XQ  ■ 40  X,  Eq.  (5) 

gives  8cr  ■ 0.72  X for  the  single  intersection  point. 

The  behavior  of  the  energy  distribution  with  film  thickness  d is 

studied  in  Fig.  5 for  the  case  that  E^  is  very  small  (E^  « E < E for 

o o — cr 

any  8).  Note,  from  (9),  that  for  8 -*■  0,  hence  y -*■  0,  A_  ■+■  Acr/2  ■ Aq  and 

E -*■  0.25  E . Since  8 ■ 1.15  X,  the  chosen  values  of  8 straddle  8 

— cr  cr  cr 

For  8 * 2.0  X and  8 * 1*2  X there  are  no  intersections,  and  for  8 ■ 1.0, 

0.9,  0.8,  and  0.5  X there  are  two  Intersections  apiece,  which  separate 

in  energy  as  8 gets  smaller.  Clearly,  the  no-intersection  cases,  8 * 2.0 
and  8 ■ 1.2  X,  correspond  to  runaway  of  the  entire  distribution,  exhibited 
by  a smooth,  rapid  increase  of  E with  d.  What  is  striking  about  Fig.  5 
is  the  demonstration  that  even  with  two  intersections  the  distribution 
appears  to  be  running  away,  namely,  for  3 ■ 1.0  X,  albeit  not  quite  as 
smoothly  as  for  8 > 8cr*  As  8 gets  smaller,  so  that  E+/E_  gets  larger 
(see  Table  I)  the  mean  energy  tends  more  to  'settle  down',  though  it  may 
still  Increase  slowly,  as  is  seen  comparing  the  8 ■ 1.0,  0.9,  and  0.8  X 
curves  in  Fig.  5.  Thus,  the  8 " 0.8  X curve  is  quantitatively  quite 
different  from  the  complete  run-away  curves,  8 '■  2.0  and  1.2  X.  The 


4000  8000 

FILM  THICKNESS  , 


Fig.  5:  The  average  energy  of  Che  distribution  as  a function  of  filn 

thickness  for  various  values  of  6*  In  all  cases,  1 - 25  X, 
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8 ■ 0.9  and  0.8  X curves  are  readily  Interpreted  in  terms  of  the  runaway, 
in  energy,  of  a small  numbers  of  electrons,  the  remainder  of  the  dis- 
tribution behaving  quite  stably,  as  illustrated  in  the  histogram  of 
Fig.  C.  Effectively,  we  have  in  these  situations  a bi-modal  distribution, 
one  containing  most  of  the  electrons  being  stable,  and  the  other  consis- 
ting of  an  unstable  energy  tall.  The  few  electrons  in  the  tail  are  those 
that  managed  to  get  past  E+.  Once  an  electron  is  well  past  E+,  the  like- 
lihood of  its  ever  returning  below  E+  becomes  vanishingly  small  — that  is, 
a true  runaway  has  occurred.  However,  as  8 gets  smaller,  E+/E  gets 
very  large,  as  seen  in  Table  1.  It  then  becomes  Increasingly  unlikely 
that  any  electron  will  get  up  to  E+  starting  from  Eq  near  0.  Thus,  for 
8 - 0.5  X in  Fig.  5 not  a single  electron  'escaped'  from  the  distribution, 
past  E+  ■ 17  Ecr  ■ 37  eV,  and  hence  the  distribution  appears  to  be  'totally 
stable,'  i.e.,  Eflv  is  quite  flat  with  increasing  d.  It  should  be  pointed 
out,  however,  that  although  this  'total  stability'  may  be  'practically' 
true,  it  cannot  be  rigorously  true.  For  a fixed,  small  8,  if  a large 
enough  number  of  electrons  are  tested,  i.e.,  if  the  batch  number  N is 
taken  large  enough,  there  will  always  be  a few  electrons  that  will  run 
away  in  energy.  This  was  simply  not  seen  for  6 - 0.5  X in  Fig.  5 because 
N - 200  was  apparently  not  large  enough.  Even  a single  electron  running 
away  for  fixed  N will  prevent  Eftv  from  being  completely  flat  with  increasing 
d,  as  seen  in  the  8 ■ 0.8  and  0.9  X curves.  The  case  8 ■ 1.0  X in  Fig.  5 
is  a transition  case,  in  which  E+  is  close  enough  to  E that  large  numbers 

of  electrons  are  running  away  in  energy.  (From  Table  I,  E+  - 2.6  Ecr  - 
5.7  eV  and  E_  ■ 0.53  Ecr  ■ 1.2  eV) . Clearly,  as  8 increases  from  1.0  X, 
through  8cr  ■ 1.15,  to  1.2  X,  the  corresponding  Eay  vs.  d plots  change 
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TABLE  I 


6(&) 

y ■ P/ecr 

+ cr 

EjE 
+ cr 

\ /X 
- cr 

E /E 
- cr 

E+/E_ 

0.5 

0.43 

4.1 

17.0 

0.57 

0.33 

51 

0.8 

0.70 

2.2 

4.8 

0.65 

0.42 

11.4 

0.9 

0.78 

1.9 

3.6 

0.68 

0.46 

7.8 

1.0 

0.87 

1.6 

2.6 

0.74 

0.55 

4.7 

1.2 

1.04 

— 

... 

... 

— 



2.0 

1.74 

... 

... 

... 

... 

... 

p - 1.15  1.  A - 50  i,  E - 2.2  eV 
Kcr  * cr  ’ cr 


smoothly.  In  this  smooth  transition  regime,  a substnatial  fraction  of  the 
total  number  of  injected  electrons  runs  away  in  energy. 

From  the  arguments  presented  above  it  is  apparent  that  whether  or  not  elec- 
trons are  likely  to  run  away  in  energy  can  depend  strongly  on  Eq,  their  energy 
at  injection.  This  is  seen  strikingly  in  Fig.  7,  in  which  E is  plotted 

® V 

vs • d for  F - 106V/cm,  e . - 0.1  eV,  k ■ 0.87  and  X » 40  X.  There  are 

pn  o 

four  curves:  EQ  - 0.01  eV  and  8 “ 0.5  and  1.0  X,  respectively;  EQ  ■ 10  eV 

and  8 ■ 0.5  and  1.0  X,  respectively.  Note  that  8fir  ■ 0.72  X and  Ecr  - 
5.63  eV  for  the  degenerate,  single  intersection  point.  Since  8 ■ 1.0  X > 8cr» 
the  two  8 ■ 1 X curves  completely  run  away,  as  expected.  On  the  other  hand, 
E_j_/E_  is  so  large  for  8 ■ 0.5  X,  and  EQ  - 0.01  eV  that  not  a single  electron 
in  the  batch  of  200  runs  away,  and  the  average  energy  saturates  with  increas- 
ing film  thickness,  as  with  the  8 ■ 0.5  X curve  of  Fig.  5.  The  more  inter- 
esting situations  are  those  for  which  EQ  ■ 10  eV  and  8 < 8cr*  Thus,  for 
8 - 0.5  X,  y - 8/8cr  - 0.693  so  that  E+  - 28  eV  and  E_  - 2.34  eV.  Although 
Eq  - 10  eV  is  a fair  amount  smaller  than  E+,  it  is  close  enough  that  a few 
electrons  are  able  to  get  past  E+  and  run  away,  while  all  the  others  settle 
down  to  a much  lower  average  energy  near  E_.  Note  that  the  8 ■ 0.5  X curve 
first  'cools  down'  and  then  rises  more  slowly.  The  cooling-down  phase  is 
dominated  by  the  main  part  of  the  distribution,  which  is  settling  down 
towards  E_  ■ 2.34  eV  from  the  high  injection  energy,  Eq  ■ 10  eV.  The  effect 
of  the  few  electrons  that  are  already  running  away  is  burled  under  this 
bulk  cooling  effect.  However,  once  the  bulk  of  the  distribution  has  settled 
down,  in  about  5000  X,  the  runaway  electrons  manifest  themselves  in  s slow, 
steady  increase  in  Eay,  as  seen  in  Fig.  7 and  in  Fig.  5,  for  the  curves 

labelled  8 - 0.8  and  0.S  X.  If  E were  larger,  correspondingly  more  electrons 

o 

would  run  away.  For  » E+  the  entire  distribution  would  run  away.  On 


KINETIC  ENERGY  OF  ELECTRON  , E («V) 

Typical  histogram  for  quasi-stable  distribution, shotring  its  bl-aodal  character 

The  parameter  values  are  k ■ 0.87,  F - 106  V/cm,  e - 0.1  eV,  X - 25  X, 
o ph  o 

6 - 0.5  X,  and  Eq  ■ 10  eV.  The  film  thickness  is  4000  X.  The  batch  siae  is 

200  particles . Note  the  three  runaway  particles  at  energies  above  20  eV. 
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the  other  hand,  6 ■ 1.0  % corresponds  to  y > 1 and  there  are  no  Inter- 


sections of  the  E (X)  and  X(E)  curves.  In  this  case  runaway  of  the 

"V 1 88 

entire  distribution  Is  expected,  independent  of  EQ. 

The  existence  of  bi-modal  distributions  leads  to  the  notion  of  local 


stability  around  a point  of  intersection  between  the  X(E)  and  E (X) 


curves.  If  the  slope  of  X(E)  is  greater  than  that  of  E (X)  at  the 

8V  9 88 

point  of  intersection,  E^,  the  distribution  is  stable  about  this  point 
in  the  following  sense:  An  electron  having  energy  slightly  greater  than 


E.  has  a mfp  that  is  shorter  than  the  one  given  by  E (X)  evaluated 

X — — flV  1 88 

at  that  particular  electron  energy.  Therefore,  the  electron  tends  to 

1 I 

cool  down  toward  E. , where  the  mfp  is  that  given  by  E (X)  evaluated 

X flV } 88 

at  E^.  Conversely,  an  electron  with  energy  less  than  E^  has  a mfp  longer 

I 

than  that  given  by  E (X) , and  it  tends  to  warm  up  to  the  point  where 

8V  ) 88 

its  mfp  is  equal  to  that  given  by  E , namely,  E . . Hence,  E.  is  a 

BY } 88  X X 

stable  point  in  energy  where  electrons  with  energies  differing  slightly 


from  E^  tend  to  evolve  toward  E^.  Note  that  X - constant  is  a special 
case  where  the  slope  of  X(E)  is  infinite.  (We  consider  negative  slopes 


of  X(E)  to  be  even  larger). 

When  the  slope  of  X(E)  is  less  than  that  of  E (X)  at  E . , insta- 

av  y 88  X 

bllity  results.  An  electron  having  an  energy  slightly  above  has  a 
mfp  which  is  greater  than  that  given  by  E __  (X)  evaluated  at  the 

8V  9 88 

electron  energy.  The  electron,  therefore,  tends  to  warm  up  to  the  energy 

E (X)  corresponding  to  its  mfp.  As  it  does  so,  its  mfp  increases 

still  further,  leading  to  a runaway  in  energy.  For  an  electron  energy 

below  E. , the  mfp  is  less  than  that  given  by  E (X)  evaluated  at 

x a V 9 88 

the  electron's  energy,  so  that  the  electron  tends  to  cool  down  toward  an 
energy  below  E.  where  X(E)  » E (X).  In  the  specific  case  where 

X 8V9SB 
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X ■ Xq  + 8 (E/eph)  and  Eav,s8  ^ k (FX)2/eph,  electron®  with  energies 


less  than  tend  to  cool  down  toward  E_ • If,  by  chance,  an  electron's  1 j 

energy  rises  above  E+,  runaway  is  strongly  favored.  Thus,  bi-modal  dis- 
tributions arise , with  most  of  the  electrons  near  the  lower  point  of  Inter- 
section, E_ , and  with  a few,  having  escaped  beyond  E^,  running  away. 

While  the  notion  of  local  stability  gives  a qualitative  feeling  for 
the  generation  of  bl-modal  distributions,  their  detailed  shape.  In  energy, 
strongly  depends  on  the  relative  slopes  of  the  curves  and  the  proximity 
of  the  points  of  Intersection. 

IV.  Conclusions 

The  conclusions  that  follow  from  I and  from  the  present  paper  are 
as  follows: 

a)  For  a constant  mfp,  X,  any  degree  of  forward  scattering,  short 

of  pure  forward  scattering,  results  in  a stable,  steady-state  energy 

distribution.  For  an  energy-dependent  mfp,  X-X(E),  our  results  are  most 

conveniently  expressed  graphically  In  the  E (ordinate)  , X (abscissa) 

plane  In  terms  of  the  intersections  of  the  two  curves  E ■ E (X) I 

av,ssv  7 , 

I X“const 

and  X - X (E) . 

b)  If  there  is  a single  Intersection  point  (which  is  not  a degeneracy), 

then  a stable  distribution  results,  with  E given  by  the  E co-ordinate 

of  the  intersection  point. 

c)  If  there  are  no  intersection  points  the  distribution  runs  away 
In  energy. 

d)  If  there  are  two  Intersection  points  the  distribution  is,  at  best, 
quasi-stable.  That  Is,  if  the  two  intersection  points  are  far  enough  apart 


in  energy,  E+  » E^,  and  Eq  «,  E+,  the  bulk  of  the  distribution  will 
be  stable,  but  there  nay  be  a few  runaway  electrons.  The  closer  Eq  gets 
to  E+  the  more  runaway  electrons  there  will  be.  In  limiting  cases,  with 
both  E_  <«  E+  and  Ec  «<  E+  the  distribution  may  be  stable,  practically 
speaking;  that  is,  the  fraction  of  the  total  number  of  particles  running 
away  may  be  too  miniscule  to  have  any  experimental  significance. 

e)  More  generally,  an  Intersection  point  at  which  A«*A(E)  has  the 
steeper  slope  [A  ■ constant  has  infinite  slope;  negative  slope  is  to  be 
considered  even  steeper]  promotes  local  stability  in  its  vicinity.  Con- 
versely, an  intersection  point  at  which  E (X)  has  the  steeper  slope 

&V  9 88 

promotes  local  runaway  in  its  vicinity.  However,  the  proximity  of  the  two 
intersection  points  is  also  significant  in  determining  the  entire  distribution. 

Finally,  we  should  like  to  stress  that  virtually  all  of  our  studies 
have  been  based  on  relatively  simple  forms  of  angular  scattering.  We  are 
simply  not  able  to  say  whether  more  complicated  scattering  processes  can 
lead  to  additional  interesting  results. 
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